Source-free electromagnetism's canonical fields 
reveal the free-photon Schrodinger equation 

Steven Kenneth Kauffmann 
American Physical Society Senior Life Member 

43 Bedok Road 
#01-11 

Country Park Condominium 

Singapore 469564 
Tel & FAX: +65 6243 6334 
Handphone: +65 9370 6583 

and 

Unit 802, Reflection on the Sea 
120 Marine Parade 
Coolangatta QLD 4225 
Australia 
Tel/FAX: +61 7 5536 7235 
Mobile: +61 4 0567 9058 

Email: SKKauffmann@gmail.com 
Abstract 

Classical equations of motion that are first-order in time and conserve energy can only be quantized 
after their variables have been transformed to canonical ones, i.e., variables in which the energy is the 
system's Hamiltonian. The source-free version of Maxwell's equations is purely dynamical, first-order in 
time and has a well-defined nonnegative conserved field energy, but is decidedly noncanonical. That should 
long ago have made source-free Maxwell equation canonical Hamiltonization a research priority, and 
afterward, standard textbook fare, but textbooks seem unaware of the issue. The opposite parities of the 
electric and magnetic fields and consequent curl operations that typify Maxwell's equations are especially 
at odds with their being canonical fields. Transformation of the magnetic field into the transverse part of 
the vector potential helps but is not sufficient; further simple nonnegative symmetric integral transforms, 
which commute with all differential operators, are needed for both fields; such transforms also supplant 
the curls in the equations of motion. The canonical replacements of the source-free electromagnetic fields 
remain transverse-vector fields, but are more diffuse than their predecessors, albeit less diffuse than the 
transverse vector potential. Combined as the real and imaginary parts of a complex field, the canonical 
fields prove to be the transverse-vector wave function of a time-dependent Schrodinger equation whose 
Hamiltonian operator is the quantization of the free photon's square-root relativistic energy. Thus proper 
quantization of the source-free Maxwell equations is identical to second quantization of free photons that 
have normal square-root energy. There is no physical reason why first and second quantization of any 
relativistic free particle ought not to proceed in precise parallel, utilizing the square-root Hamiltonian 
operator. This natural procedure leaves no role for the completely artificial Klein-Gordon and Dirac 
equations, as accords with their grossly unphysical properties. 



Introduction 

Notwithstanding the approximately century and a half which has passed since both the development of 
Hamiltonian classical dynamics and also the final codification of the laws which govern electromagnetic fields 
in configuration space, canonical Hamiltonian formulation of the purely dynamical source-free instance of 
electromagnetic field theory in configuration space unaccountably still lingers as essentially terra incognita. 
This is all the stranger insofar as proper canonical Hamiltonization of a classical dynamical system is a 
fundamental prerequisite to its rigorous unambiguous quantization, whether by the Hamiltonian phase-space 
path integral [1], or by the self-consistent (slight) extension of Dirac's canonical commutation rule [2], 

A crucial aspect of correct canonical Hamiltonization of the source-free instance of Maxwell's field equa- 
tions in configuration space is that it is subject not only to those equations themselves, but also to the 
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fact that the conserved energy content of the source-free (transverse) dynamical electromagnetic fields is 
a very specific known functional of those fields. This specific conserved field-energy functional necessarily 
becomes the field system's Hamiltonian functional when its fields have been transformed to canonical ones, 
i.e., we shall have transformed the fields to canonical ones when the field equations expressed in terms of 
the transformed fields exactly match the Hamiltonian equations of motion for the transformed fields which 
result from taking the Hamiltonian to be the conserved field-energy functional as expressed in terms of the 
transformed fields. Therefore the fact that the conserved field-energy functional is known a priori places a 
strong constraint on just what the canonical fields can be. 

In the source-free case, Maxwell's field equations in configuration space are strictly linear and homoge- 
neous, so it is not surprising that we can restrict ourselves to purely linear transformations of the dynamical 
magnetic and transverse electric fields in the search for properly canonical fields. (The longitudinal elec- 
tric field is never dynamical in character, and it vanishes identically in the source-free case.) It as well 
turns out to be unnecessary to mix the magnetic field with the transverse electric field in the course of this 
search. But the axial/polar dichotomy between the magnetic field and the transverse electric field turns out 
to be incompatible with fields that are properly canonical; this failing is readily rectified by a one-to-one 
orthogonal linear mapping of the axial- vector magnetic field onto a transverse polar- vector counterpart. The 
resulting transformed field equations, however, still fail to match the Hamiltonian equations of motion which 
follow from the transformed conserved field-energy functional, but a further nonncgativc symmetric linear 
transformation of all fields, which commutes with all partial derivative operators and doesn't affect the field 
equations, repairs this disparity. This last transformation is decidedly nonlocal however: the still transverse 
canonical fields turn out to be diffused relative to the transverse electric and magnetic fields — albeit to 
a lesser extent than the transverse vector potential is likewise diffused relative to the magnetic field (the 
Aharanov-Bohm effect is a quite dramatic manifestation of that particular relative diffusion). 

The transverse canonical vector fields are not only diffused relative to the transverse dynamical electric 
and magnetic fields; they also have a different dimensionality, namely that of the square root of action density 
rather than that of the square root of energy density, which is the dimensionality of electric and magnetic 
fields. If the two transverse canonical vector fields are divided by (271)5 and then combined as the real 
and imaginary parts of a complex transverse-vector field, it is seen that this complex-valued field satisfies 
the natural time-dependent Schrddinger equation in configuration representation whose Hamiltonian oper- 
ator is hc(— V 2 )2, namely |cp|, the Hamiltonian operator for a first-quantized free solitary ultrarelativistic 
zero-mass particle that is the zero-mass limit of the natural correspondence-principle-mandated square-root 
relativistic free-particle Hamiltonian (m 2 c 4 + cp | 2 ) 5 . Moreover, the dimensionality of the complex- valued 
transverse- vector field in this Schrddinger equation is that of the square root of probability density, which is 
appropriate for a Schrddinger equation wave function. In other words, the two transverse canonical vector 
fields of source- free electromagnetic field theory comprise the natural basis for the complex- valued transverse- 
vector wave-function description of the first- quantized free solitary photon. Source-free electromagnetic field 
theory is in this way precisely the time-dependent Schrddinger-equation description of the first-quantized 
free solitary photon with its natural square-root Hamiltonian operator, but that stark fact simply does not 
become technically manifest until this field theory has been properly canonically Hamiltonizedl 

It is delightfully amazing that James Clerk Maxwell, building atop the foundation laid by Michael Fara- 
day's experimental results, effectively discovered the Schrddinger equation, indeed for a tricky transversely- 
polarized spin 1 ultrarelativistic particle, very long before Erwin Schrddinger's own nonrelativistic quantum 
insights. Remarkably, Maxwell could accomplish this feat with no knowledge whatsoever of the quantum 
of action — whose discovery by Max Planck still lay well in the future — because of the intriguing happen- 
stance that in configuration representation the photon's zero-mass property "releases" a factor of h from its 
Hamiltonian operator that neatly cancels out the factor of U which is a fixture of the left-hand side of the 
time-dependent Schrddinger equation. Thus does ultrarelativistic quantum mechanics chameleon-like meta- 
morphose into "classical field theory"! The technical details of this connection can not, of course, be fully 
laid bare until Maxwell's effective formulation of the theory in terms of the axial magnetic and transverse- 
polar electric fields has been properly canonically Hamiltonized — the two corresponding transverse canonical 
fields have the same (not the opposite) parity, are relatively somewhat more diffused than their antecedent 
magnetic and electric fields, albeit less so than is the case for the transverse vector potential relative to its 
antecedent magnetic field, and are joined together in a single complex- valued transverse- vector wave function 
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that describes the first-quantized free solitary photon. 

In the world of quantum theory, canonical Hamiltonization becomes merely the required prelude to 
quantization; this of course also applies to the quantization of source-free electromagnetic field theory. After 
completion of canonical Hamiltonization, canonical quantization parlays the Poisson bracket relations of the 
fundamental canonical dynamical variables into commutation relations. In canonically Hamiltonized source- 
free electromagnetic field theory, the fundamental canonical dynamical variables are the vector components 
of the two transverse-vector canonical fields themselves. Thus it is the components of these transverse- 
vector canonical fields themselves that quantization promotes into noncommuting Hcrmitian operators; these 
operators inherit from their simple mutual Poisson bracket relations equally simple mutual commutation 
relations. The vector components of the free solitary-photon wave functions, being straightforward complex 
linear combinations of the corresponding vector components of the two real transverse- vector canonical fields, 
become non-Hermitian noncommuting operators such that there are as well simple commutation relations 
between the vector components of the quantized wave function and those of its Hermitian conjugate(the 
transverse character of these canonical vector fields and complex vector wave functions subtracts an annoying 
longitudinal projection-operator term from their formal configuration-representation commutators, which 
creates an ugly and potentially confusing distraction that is merely technical in nature). The Hamiltonian 
functional also becomes quantized via its bilinear dependence on the real canonical fields (or, as well, via its 
linear dependence on both the complex photon wave function and its complex conjugate), and is a Hcrmitian 
operator. As is the normal case in quantum theory, the Heisenberg picture with respect to this Hamiltonian 
operator manifests the same dynamical equations for the quantized field operators as they obeyed prior to 
their quantization. Thus the quantized wave function operator still satisfies the very same Schrddinger 
equation that it satisfied prior to its quantization, which can thus properly be termed second quantization. 
This Schrodinger equation in its second-quantized form obviously still features the very same m — > limit of 
the relativistic square-root Hamiltonian operator (m 2 c 4 + |cp| 2 ) 2 (which comes to fic(— V 2 )^ in configuration 
representation) that it featured prior to that quantization! Finally, the simple commutation relation between 
the quantized non-Hermitian wave function operator and its Hermitian conjugate is such that the quantized 
wave function operator is intcrprctable as the annihilation operator for free photon states in the underlying 
Hilbert space (commonly called Fock space), while its Hermitian conjugate is likewise interpretable as the 
creation operator for such free photon states. 

This second- quantized theory of arbitrarily many free photons, or, equivalently , quantized theory of 
source-free electromagnetism, is not, of course, the end of the physics story. We have so far postponed 
consideration of the coupling of charged matter to the dynamical transverse part of the electromagnetic 
fields, i.e., to photons. In the rather global technical language of the Maxwell equations, this coupling (in 
fact an inhomogeneous driving term) can be mathematically abstracted in terms of the transverse part of the 
current-density vector field. From a microscopic charged-particle perspective, however, relativistic coupling 
to electromagnetism occurs via the four- vector potential A 11 . The part of which is related to the dynamical 
transverse electromagnetic fields is, naturally enough, the transverse part of the vector potential A, which we 
denote as At- Indeed the precise linear relationship of At to the complex- valued transverse- vector photon 
wave function and its complex conjugate can be readily traced. The remaining A and longitudinal part of 
A are related 1) to the nondynamical longitudinal part of the electric field, which is a mere homogeneous 
functional of the global charge density, and 2) to a nonphysical "gauge degree of freedom" . Thus they bear 
no direct relation to the photon wave function. Because the photon wave function and its complex conjugate 
become operators upon second quantization, At becomes a linearly related operator as well, in fact a 
Hermitian one that can either create or annihilate a free photon state in the underlying Fock space. This 
detailed tying of At to the free-photon annihilation and creation operators is the technical way in which the 
direct interaction between photons and charged particles finds effect in interacting second-quantized theory. 
Such a theory must of course as well includes the above-described second-quantized Hamiltonian functional 
operator of source- free electromagnetism, i.e., the Hamiltonian operator for the for an arbitrary number of 
nonintcracting free photons. 

What about ^4° and the longitudinal part of A, which we write as A^? These are also coupled to charged 
particles, but have no direct relation to dynamical transverse electromagnetism or photons. They encompass 
the nondynamical physics of E^, the longitudinal part of the electric field, which can be deduced from 
Maxwell's equations to be entirely the creature of the global charge density, and therefore is properly termed 
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coulombic, and the "nonphysics" of the choice of gauge function. Given the global Hamiltonian functional for 
the charged particles, specifically including in particular their formal interaction with electromagnetism via 
arbitrary , we readily obtain the global charge density p as the functional derivative of that Hamiltonian 
functional with respect to A . The longitudinal part of the electric field, namely E^, is then obtained as a 
closed- form linear homogeneous functional of p, but this is not sufficient to pin down both A and Al, which, 
additionally, calls for a choice of gauge. The simplest choice is the Coulomb gauge, which rather brusquely 
puts Al to zero, oblivious to even the slightest pretense of special- relativistic finesse. Thereupon A also 
becomes a closed-form linear homogeneous functional of p via the static Coulomb kernel. Aside from the 
blatantly nonrelativistic character of a potential A that arises from the static Coulomb kernel, two other 
subtleties present themselves: 1) since A couples to the very same particles that give rise to the global p, 
one must compensate for double counting by the usual expedient of halving A relative to its naive value, 
and 2) if it should turn out that p itself still retains a formal dependence on A , one will have only obtained 
an implicit equation for ^4°, one that requires solution (probably only an approximate one via successive 
iteration will be feasible). 

A relativistically more plausible gauge [3] is rooted in the Lorentz condition (cV • Al + A ) =0 conjoined 
to the requirement that A be linear and homogeneous in p, but retarded by (|r — r'|/c) This can be achieved 
because the Lorentz condition imposes on A the requirement that it be related to p by the c-speed compatible 
second-order partial differential equation (A°/c 2 — V 2 A°) = p. The upshot is that ^4° is a closed- form integral 
that is linear and homogeneous in p via the static Coulomb kernel, but for p evaluated at the retarded time 
it— |r— r'|/c). With A in hand, Al can be determined from ^4° via the Lorentz condition. In the special case 
that p = 0, this "retarded Lorentz gauge" is identical to the Coulomb gauge. The above caveats concerning 
halving ^4° relative to its naive value and the possibility of obtaining only an implicit equation for A still 
apply. 

The following sections provide the technical/mathematical details that explicitly underpin and demon- 
strate what is described and discussed in the foregoing paragraphs. 



Longitudinal nondynamical and transverse dynamical electomagnetic fields 

The four Maxwell equations for the electromagnetic field (E, B) with four-current source (p,j/c) are com- 
prised of Coulomb's law, 

V • E = p, (la) 

Faraday's law, 

VxE=-B/c, (16) 

Gauss' law, 

VB = 0, (lc) 

and Maxwell's law, 

VxB = (j + E)/c. (Id) 

The Coulomb and Gauss laws do not involve time derivatives of the electromagnetic fields, i.e., they are 
nondynamical in character. This raises the possibility that some part of the electromagnetic field (E, B) may 
itself be of nondynamical character, i.e. determinable without reference to any initial conditions. This in 
fact turns out to be the case for the longitudinal part of the electric field, which it is therefore very useful to 
detach from the rest of the electromagnetic field. The ability to separate a vector field into its longitudinal 
and transverse parts in a linear and unique fashion is in fact extremely useful throughout electromagnetic 
theory, so we turn first to a discussion of how that is carried out. 

We shall now indicate that any vector field F(r) which is continuously differentiable, and for which 

V • F(r)|/|r| is intgrable over all space, has a unique decomposition achieved by linear operation into the 
sum of its longitudinal part Fz,(r) with its transverse part Fr(r), where V • F^ = V • F, VxFi = and 

V • F T = 0. 

Since we require that V x F^ = 0, it will need to be the case that F^(r) = — V5(r), where S(r) is a 
scalar field. Since we also require that V • F^ = V • F, we must have that — V 2 5 = V • F, whose general 
solution is, 

5(r) = c + k .r+[(-V 2 )- 1 (V-F)](r), 
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where Co and ko are arbitrary constants, and the natural notation (—V 2 ) 1 denotes the integral operator 
whose configuration- representation kernel (r| (— V 2 ) -1 |r') is, 

(r|(-V 2 )- 1 |r') = (4 7 r|r-r'|)- 1 , 

because, as is well-known, 

(-V 2 )(47r|r - r'l)- 1 = <5( 3 )(r - r') = (r|r'}. 

Since Fl = — VS, we obtain Fl = ko — V[(— V 2 ) _1 (V • F)], where ko is an arbitrary constant. Because 
we require the decomposition of F into the sum of Fl with Ft to be achieved by linear operation, k must 
vanish identically, and therefore, 

F L = -V[(-V 2 )- 1 (V-F)], 

or 

F L (r) = -V r /(4rr|r - r'l)" 1 ^ • F(r')) dV. 

We now note from the above two equations the key fact that Fl is entirely determined by V • F. Since, of 
course, Ft = F — Fl, the fact that V • Fl = V • F implies that V • Ft = 0, as is required for Ft- 

We shall now systematically apply this linear decomposition into its longitudinal and transverse parts to 
each vector field in each of the four Maxwell equations of Eqs. (1). We consequently obtain, 

V • E L = p, (2a) 

V x E T = -B T /c, (2b) 

B L = 0, (2c) 

Bl = 0, (2d) 

V x B T = (jt + E T )/c, (2e) 

and, 

jl + El = 0. (2/) 

From Eq. (2a) and the representation for the longitudinal part of any vector field obtained above, we obtain 
El in closed form, 

E L (r, t) = -V r [(-V 2 )-V](r) = -V r /(47r|r - r'|)- V(r', t) dV. (3a) 

Eq. (2d) makes Eq. (2c) redundant. In consequence of Eq. (2d), B = Bt, which we shall simply always bear 
in mind. That enables us to drop all references to B T - Therefore Eq. (2b) can be written, 

V x E T = -B/c, (36) 

and Eq. (2e) can likewise be written, 

V x B = (jT + E T )/c. (3c) 

With regard to Eq. (2f), we recall that the longitudinal part of a vector field is completely determined by its 
divergence. Therefore it is sufficient to simply work with the divergence of Eq. (2f), from which we deduce 
that V • j + V • El = 0. Now there is nothing more to be learned about El, as it is given in closed form in 
terms of p in Eq. (3a). Therefore we use Eq. (2a) to eliminate its presence, and thus obtain the celebrated 
constraint due to charge/current conservation: 

V-j + p = 0. (3d) 

Eqs. (3), plus the fact that B is purely transverse, completely replace the Maxwell equation system of Eqs. (1), 
and are clearly far more informative than that system. In particular, it is crystal-clear from Eq. (3a) that 
El is a completely nondynamical variable which is entirely independent of the choice of initial conditions. It 
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is only the transverse fields Et and B that are actual dynamical variables, and only these can legitimately 
be incorporated into a standard dynamical framework! 

Source-free electromagnetic theory 

We shall now do away with the charged matter sources, namely put p and j to zero, and thereby deal with 
purely self-sustaining radiation. Eq. (3a) thereupon becomes E^ = 0, and thus can be dropped entirely along 
with Eq. (3d), which reduces to the trivial identity — 0. We are left with the following two dynamical 
equations, which involve only the two purely transverse fields Et and B, 

B = -cV x E T , (4a) 

and, 

Et = cV x B. (46) 

In addition to these equations of motion, source-free electromagnetism has a very well-known conserved 
nonnegative field-energy functional [4,5], which is given by, 

£[E T ,B] = 1/ [|E T (r,t)| 2 + |B(r,t)| 2 ] d 3 r. (4c) 

We can readily calculate the time rate of change of E[J$t, B] by applying the two field equations of motion, 
namely Eqs. (4a) and (4b), which yields, 

dE/dt = cf [(V x B) • E T - B • (V x E T )] d 3 r (Ad) 

Now it is an identity that, 

V • (B x E T ) = (V x B) • E T - B • (V x E T ), 

which implies that, 

dE/dt = c j V • (B x E T ) d 3 r, (Ae) 

and the integral over all space of such a divergence will, of course, vanish under normal circumstances. Thus 
the nonnegative field-energy functional £?[Et,B] is indeed conserved, 

dE/dt = 0. (4/) 

A useful corollary of this demonstration is that under normal circumstances, 

/(V x F) • G d 3 r = J F • (V x G) d 3 r. (5) 

Now the correct conserved energy of a dynamical system becomes the system's Hamiltonian whenever that 
system is described by correct canonical variables. Here it is immediately clear that (Et, B) are not correct 
canonical fields for this dynamical electromagnetic system, because treating the system's correct energy 
functional of Eq. (4c) as the system's Hamiltonian functional results in the putative Hamiltonian equations 
of motion El = B and B = — E^, which blatantly disagree with the actual equations of motion that are 
given by Eqs. (4b) and (4a) — there is no trace of the very prominent curl operations of the equations of 
motion to be found in the much more austerely straightforward Hamiltonian equations that flow from the 
nonnegative field-energy functional. 

Like parities and additional diffuseness of the canonical fields 

We are now clearly obliged to search for correct canonical fields for this source-free electromagnetic system 
by trying out transformations of (Et, B). Since we are dealing with linear equations and a bilinear energy 
functional, it is clear that we can restrict ourselves to linear transformations. Both the form of the bilinear 
energy functional and that of the present equations of motion of Eqs. (4b) and (4a) strongly suggest that 
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we not look at any transformations that mix the electric with the magnetic field. In light of the form of 
the bilinear energy functional, however, it seems urgent to find a transformation that eliminates the curl 
operations from the equations of motion. The presence of the curl operations is entwined with the fact that 
B is an axial transverse vector while is a polar transverse vector; thus it seems important to find a way 
map B onto a polar transverse vector without harming any physically important information carried by B. 

The vector potential A is a polar vector whose tranverse part appears to carry all the physical information 
that is contained in the transverse axial vector B. This is so because B = VxA = Vx At- However, At 
itself has a different dimensionality from B, and is also considerably more diffuse, as the Aharanov-Bohm 
effect pointedly illustrates. It would be good to be able to apply a transformation to At that leaves its 
transverse polar nature intact, yet compensates for its differences from B in dimensionality and diffuseness. 
In this regard, it is interesting to note that V x B = V x (V x At) = — V 2 At, which is less diffuse than 
B, and errs in dimensionality relative to B in the opposite direction from the dimensionality error made by 
At- These considerations strongly suggest that the object we would truly like to have is (— V 2 )^ At- Since, 
as we have noted just above, V x B = — V 2 At, we have that, (— V 2 )^ At = (— V 2 )~2 V x B, i.e., we can 
conveniently express the entity we want entirely in terms of B itself, without any need to make reference 
to At- It is convenient to invent a shorthand notation for this desired transformation of B from an axial 
transverse vector to a polar transverse vector that otherwise mimics B itself just as closely as possible, 

B* = f (— V 2 )~5 V x B, (6a) 

which has the marvelous property that, 

(B*)t = B, (66) 

i.e., this linear operation is actually a conjugation when it is restricted to transverse vector fields . We can dub 
it "polar-axial conjugation" . To quell any lingering doubts as to the calculational "meat" on its somewhat 
abstract "bone", we explicitly exhibit the integral-operator kernel of (— V 2 )~2 in configuration represention, 

<r|(-V 2 )^ |r'} = (2^)- 3 | |k|- 1 e' k -( r - , )rf 3 k= 27r2|r 1 _ r/|2 . (6c) 

It is, moreover, very satisfying to note that the physically key transverse part of the vector potential is neatly 
related to this conjugate of B, 

A r = (-V 2 )-5B*. (6d) 

Finally, just as a conjugation operation ought to be, polar-axial conjugation is orthogonal in the natural 
Hilbert space of transverse vector fields . Thus the conserved nonnegative field-energy functional E[Et, B] of 
Eq. (4c) is invariant under this conjugation. In particular, from the corollary given by Eq. (5), the symmetric 
nature of the linear operator (— V 2 )~2 ; the fact that it commutes with differential operators and the fact 
that its square equals (— V 2 ) -1 , it follows that, 

/ B* • B*d 3 r = / B • V x [V x (-V 2 )- X B] d 3 r = / B • B d 3 r. (6e) 

Although the polar-axial conjugation of B leaves the conserved nonnegative field-energy functional 
£7[Et,B] of Eq. (4c) form-invariant, it does abolish the curl operations from the equations of motion. 
Thus Eq. (4b) can immediately be rewritten in terms of B* as, 

E T = c(-V 2 )5B*. (7a) 

Writing Eq. (4a) in terms of B* involves slightly more work, in that the curl operator must first be applied 
to both sides before the translation in terms of B* can be made, 

B* = -c(-V 2 )*E T . (76) 

The operator (— V 2 )^ is most easily handled in Fourier transform, as it is distribution- valued (has locally 
singular behavior) in configuration representation. It's kernel there is given by, 

(r|(-V 2 )^|0 = (2 7 r)- 3 / |k|e^-Vk=Hm j^j^j^ . (7c) 
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Having knocked the curl operations out of the equations of motion, Eqs. (7a) and (7b), we are much 
closer to our goal of achieving canonical fields, but are not quite there yet. However, what still needs to 
be done is now relatively easy to light on. The operator factor with dimensionality of frequency that now 
uniformly emerges in both equations of motion, namely c(— V 2 )^ must also be persuaded to show its face in 
the transformed field-energy functional: our polar-axial conjugation transformation changed the equations 
of motion and left the field-energy functional form-invariant, but our next transformation must do precisely 
the opposite. Multiplying both of and B by a power of the operator [c 2 (— V 2 )] doesn't change the form 
of Eqs. (7a) and (7b) at all, but it causes that operator to the negative of twice that power to appear in the 
field-energy functional. To have canonical consistency between the equations of motion and the field-energy 
functional, the power of the operator [c 2 (— V 2 )] that must appear explicitly in the field-energy functional is 
one half. That requires the canonical fields to be equal to the present E T and B* fields times this operator 
to the power of minus one quarter. 

Our canonical fields are now "in the bag" . They are explicitly, 

* = [c 2 (-V 2 )]-JE T , n= [c 2 (-V 2 )]-iB*. (8a) 

The equations of motion don't change their form from that of Eqs. (7a) and (7b), notwithstanding our having 
changed to the canonical fields <I> and n, because the common operator factor of [c 2 (— V 2 )] 3 simply factors 
through and out of those equations. Thus the equations for <j> and II are, 

* = c(-V 2 )5n, 11= -c(-V 2 )**. (86) 

When the field-energy functional is written in terms of the canonical fields <I? and n, it does change its form; 
it now becomes a Hamiltonian that is consistent with the equations of motion of Eq. (8b). This Hamiltonian 

is, 

#[*,II] = if * (c(-V 2 )^*) + 11 (c(-V 2 )^n) d 3 r. (8c) 

We noted in Eq. (6d) that A T = (-V 2 )~i B*. Now from Eq. (8a) we deduce that B* = c^(-V 2 )3lI. 
Therefore we obtain that, 

A r = c5(-V 2 )-'n (8d) 

Eq. (8a) shows that the canonical fields are somewhat more diffuse than the corresponding electric and 
magnetic fields, and while those fields have dimensionality of the square root of energy density, the canonical 
fields have dimensionality of the square root of action density. 



The complex-valued photon wave function and Schrodinger equation 

Now let us make the transition to the standard complex transverse vector field which has the dimensionality 
of the square root of probability density. This field is, 

¥ = (n-i*)/(2&)i. (9a) 

We note from Eq. (8b) that, 

(ri - i*) = c(-v 2 )i (-* - in) = -ic(-v 2 )5 (n - 

Multiplying both sides of the above result through by ih~, we deduce from Eq. (9a) that, 

ih& = 7ic(-V 2 )5*. (96) 

Eq. (9b) is a Schrodinger equation with square-root Hamiltonian operator h~c{— V 2 )^ = |cp|, which is that of 
a massless free particle. This Schrodinger equation's wave function \& has the dimensionality of the square 
root of probability density, and it is a complex transverse vector field. In other words, after its proper 
canonical Hamiltonization and field complexification, source-free electromagnetism perfectly describes the 
solitary, first-quantized free photon, i.e., it is revealed to be first-quantized photodynamics. Let us now 
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tie down the last detail of this identification by rewriting the canonical Hamiltonian functional H[&, II] of 
Eq. (8c) in terms of \I> and its complex conjugate The result is, 

H[V, *] = /** ■ (hc(-V 2 )i *) d 3 r. (9c) 

This is indeed the Hamiltonian functional that corresponds to the configuration-space Schrodinger equation 
of Eq. (9b), whose complex- valued transverse- vector wave function \I> is now completely ready for second 
quantization. Before carrying this out, let us express At, the transverse part of the vector potential, in 
terms of the complex- valued photon wave function ty. To do this we combine the result of Eq. (8d) with 
n = {h/2)i (* + **), which follows from Eq. (9a), to obtain, 

A T = (ftc/2)* (-V 2 )-^ (* + **). (9d) 

This reveals the transverse part of the vector potential to be somewhat more diffuse than the photon wave 
function, which itself is revealed by Eqs. (9a) and (8a) to be equally more diffuse than the electromagnetic 
fields. 

Second-quantized photodynamics 

If the photon wave function 'S'(r) were a complex- valued full-vector field, it would be quantized by being 
changed into a non-Hermitian operator \& (r) whose commutation relation with its Hermitian conjugate \l>t (r) 
would be given by, 

[{9{T)) i ,{&{T>)) j ]=6 ij 5<- 3 XT-T>) = (T\8 ij \i>), i,j = 1,2,3, (10a) 

which is the straightforward consequence of Hermitian quantization of the real canonical vector fields <&(r) 
and II(r) that accords with the Dirac relation of commutators to classical Poisson brackets for the quan- 
tized components of the classical phase-space vector. Here, however, we must contend with the techni- 
cal/mathematical annoyance imposed by the transverse-vector character of the real canonical fields $(r) 
and II(r), which effectively removes one third of the classical phase-space degrees of freedom that would 
have been present had 3>(r) and II(r) been real canonical full-vector fields. What is missing from the clas- 
sical full-vector field phase space is, of course, its longitudinal part, i.e., all those vector fields which can be 
written as the gradient of a scalar field. The ij components of the projection operator into this longitudinal 
part of the space of vector fields are given by the integro-differential operators (— X7 2 )~ 1 dj). This opera- 
tor is readily verified to map any vector field into the gradient of a scalar, to be a symmetric linear operator 
on the natural Hilbert space of real vector fields and to be equal to the square of itself. Therefore the ij 
components of the projection operator into the transverse part of the space of vector fields are given by 
the operators (Sij + 9,(— V ) dj). Thus to be mathematically consistent with the transverse-vector nature 
of the classical canonical and quantum operator fields, we need to replace Eq. (10a) by the commutation 
relation, 

[(*(!•))<, (§t(r')),] - <r|(<% +3 i (-V 2 r 1 3 J )|r') = 

(106) 

M (3) (r - r') - (2tt)- 3 Je ik -( r - r ')(k) i |k|- 2 (k) j d 3 k , i,j = 1, 2, 3. 

This commutation relation has the well-known interpretation that ^(r) creates a free-photon state localized 
at r, while SJ'(r) annihilates such a free-photon state, so that the underlying Hilbert space, called Fock space, 
accommodates arbitrarily large numbers of free photons [6] . The real- valued bilinear Hamiltonian functional 
H of Eq. (9c) must be correspondingly quantized to become the Hermitian Hamiltonian operator 

H[W ,4?] that governs the time evolution of this second-quantized photodynamical system of arbitrarily- 
large numbers of free-photons. An very important (and normally expected) consequence of this is that in the 
Heisenbcrg picture the non-Hermitian field obeys the selfsame Schrodinger equation of Eq. (9b) that was its 
equation of motion before it was second-quantized. In other words, the occurrence at the quantized- field level 
of the free-photon Schrodinger equation of Eq. (9b) is an unavoidable consequence of correct quantization of 
the source- free Maxwell equations. 
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At this second-quantized level, At, the transverse part of the electromagnetic potential, also becomes 
an operator, in fact a Hcrmitian one. From Eq. (9d) we see that its operator form will be given by, 

A T = (hc/2)i(-V 2 )-i (* + *t). (10c) 

It is apparent that At is the object which couples photons to charged particles. From its form as given 
by Eq. (10c) it is clear that charged particles can both absorb and emit photons via At- We now cast a 
backward glance at what has been accomplished in the foregoing sections, and then take to its completion 
the discussion we have just begun on how to theoretically set up the electromagnetic interactions that affect 
charged particles. 

Photodynamical and coulombic interactions of charged particles 

The source-free case of electromagnetism has, after its meticulous canonical Hamiltonization, which is an 
absolute necessity that has been essentially universally honored only in the breach for the past century and 
a half, effortlessly yielded up the completely natural first- and second-quantized theories of free transverse- 
vector photons, with the later governed by the bilinear-field Hamiltonian operator iJ^^], which is 
the straightforward quantized version of the c-number Hamiltonian functional iJ[\J r *,\J'] of Eq. (9c). The 
quantized-field Hamiltonian H[W , \P] yields in the Heisenberg picture an equation of motion for the quan- 
tized field \I/ that is form-identical to Eq. (9b), the Schrddinger equation for the first-quantized photon 
wave function This second- quantized form of the Schrodingcr equation for \& obviously still features the 
very same m — > limit of the relativistic square-root Hamiltonian operator (m 2 c 4 + cp | 2 ) 5 (which comes 
to hc(— V 2 )s in configuration representation) that it featured as it applied to Thus the emergence of 
this free-photon Schrodingcr equation with the just-mentioned natural square-root relativistic Hamiltonian 
operator is the unavoidable consequence of quantization of source-free electromagnetism! 

Maxwell's equations indicate that charged matter interacts with the dynamical transverse part of elec- 
tromagnetism via a global transverse current density, as is seen from Eq. (3c). Microscopic theories of 
charged particles, however, indicate that their interaction with transverse electromagnetism occurs via At, 
the transverse part of the electromagnetic potential, whose quantized form At is given directly in terms of 
the Hcrmitian part of the non-Hcrmitian quantized photon wave-function field \I/ by Eq. (10c). 

To describe charged particles and electromagnetism participating in mutual interaction, we require a 
Hamilton which consists of the sum of with a Hamiltonian that describes the charged parti- 

cles, specifically including the effects on them which an arbitrary external c-number four-vector potential 

= (A , A) produces. After writing A = A^ + At by applying the now familiar unique linear longitu- 
dinal/traverse decomposition of any vector field, we concretely replace At by the Hermitian photon field 
operator At that is given by Eq. (10c). The remaining A° and A^ have no direct relation to photons: in 
light of the basic definition = (— \7A° — Al/c) and Eq. (3a), they are determined 1) by the global charge 
density p of the system, which exerts a coulombic effect on itself, and 2) by one's choice of gauge. The sys- 
tem's global charge density p is given by the functional derivative of the charged-partiele Hamiltonian with 
respect to A . If it should happen that this p itself depends on A , then, after making the choice of gauge, 
one will still not have in hand the actual result for A , but only an implicit equation for that result, whose 
solution can probably only be successively approximated by iteration. So far as choice of gauge is concerned, 
by far the simplest is the Coulomb gauge, which, via its requirement that V • A = 0, neatly implies that Al 
vanishes, albeit this is prima facie disrespectful of special relativistic precepts. That notwithstanding, there 
does not seem to be any physical reason to eschew the Coulomb gauge, since Maxwell's equations yield the 
result of the coulombic effect to be the Eq. (3a) form of in terms of p, which is very closely related to 
the Coulomb gauge result below for A , when A is additionally required to be linear and homogeneous in 
p, namely, 

yl (r,t) = i/(47r|r-r'|)-V(r',t)dV. (11) 

An unusual factor of one half has been inserted into Eq. (11) to compensate the double counting that would 
otherwise occur because this particular A inherently interacts coulombically with the very same charge 
density p that gives rise to it. 
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If, in spite of the very close relationship of the Coulomb gauge result of Eq. (11) to the Maxwell equation 
result of Eq. (3a) for E^, a relativistically more plausible gauge is nonetheless thought to be desirable, the 
retarded Lorentz gauge of Ref. [3] would seem to be an excellent choice. Scrapping the Coulomb gauge 
requirement that = 0, but imposing the relativistically impeccable Lorentz condition (cV • + A ) = 
implies that A satisfies the second-order in time partial differential equation (A°/c 2 — \7 2 A°) = p. Imposing 
the further requirements that A be linear and homogeneous in p, and that it respond to changes in p only 
after the c-speed retardation time |r — r'|/c, turns out to yield an A which is different from that of Eq. 11 
only in that p(r',t) on its right-hand side is replaced by p(r',t — |r — r'|/c). With that retarded-Lorentz- 
gauge result for ^4° in terms of p in hand, one then determines the retarded-Lorentz-gauge A^ from the 
retarded-Lorentz-gauge A and the Lorentz condition, which explicitly yields A^ = V[(— V ) (A /c)]. For 
the case that p is time-independent, this full retarded-Lorentz-gauge result reduces to that of the Coulomb 
gauge. 

Conclusion 

It has recently been strenuously argued on the basis of the correspondence principle that the only physically 
sensible Hamiltonian operator for a solitary, relativistic first-quantized free particle of positive mass m is the 
square root operator (m 2 c 4 + |cp| 2 )a [7]. The extension of this idea to massless particles of course assigns 
them the Hamiltonian operator |cp|, which is %c{— V 2 )^ in configuration representation. One naturally 
turns to a long and firmly established theory, namely source-free electromagnetism, which is supposed to 
include a massless free particle, namely the free photon, in its ambit, to see how this Hamiltonian-operator 
cum Schrodinger equation idea fares in the context of its quantization. The foregoing work shows that it 
fares absolutely brilliantly, with the expected Schrodinger equation and its associated Hamiltonian operator 
being perfectly maintained right to the level of the quantized field (or, more precisely, quantized photon wave 
function) if one does not neglect to properly canonically Hamiltonize the source-fee Maxwell equations (whose 
electric and magnetic fields are very far from being properly canonical!) before embarking on quantization. 
In fact, this Schrodinger equation's appearance at the first- quantized level turns out to be no more than a 
simple, direct by-product of merely the proper canonical Hamiltonization of Maxwell's source-free equations! 

These results for source-free electromagnetism lend impressive support to the almost ridiculously straight- 
forward idea that the correspondence-principle-mandated relativistic free-particle square-root Hamiltonian 
operator (m 2 c 4 + |cp| 2 ) 2 and associated time-dependent Schrodinger equation is without exception the correct 
starting point for relativistic quantum theory (a glance at the section just above, and perhaps even more so 
at the latter sections of Ref. [7] , reveals a taste of the mind-boggling complexity and richness which this mere 
starting point rapidly gives way to) . Who would be prepared to for an instant contest the completely paral- 
lel assertion that the correspondence-principle-mandated nonrelativistic free-particle kinetic-energy Hamil- 
tonian operator |p| 2 /(2m) and associated time-dependent Schrodinger equation is without exception the 
correct starting point for nonrelativistic quantum theory? It may not, in this very regard, have escaped the 
reader's attention that the correspondence-principle-mandated square root relativistic Hamiltonian has the 
delicately subtle property that as c — > 00, 

[(m 2 c 4 + \cp\ 2 )i - mc 2 ] -> |p| 2 /(2m). 

For the linearized Dirac Hamiltonian a ■ pc + lime 2 , there is simply no remotely similar property. The 
very best in this regard which can be done with the Dirac Hamiltonian is to subtract away its value at 
p = 0, leaving a ■ pc, which diverges as c — > 00! Klein-Gordon theory fails altogether to present us with 
a Hamiltonian (which is the proximate cause of its quantum-theoretic downfall [7]), but it does give us 
the square of the correspondence-principle-mandated relativistic one, namely (m 2 c 4 + cp | 2 ) . Subtracting 
away its value at p = leaves cp| 2 , which also diverges as c — > 00! The relativistic free-particle square 
root Hamiltonian (m 2 c 4: + |cp| 2 ) 2 is very subtly tuned indeed! The square root itself is in fact theoretically 
completely entwined with the square roots that are archetypal of the Lorentz transformation [7] . 

In contrast to this square root, which is uniquely fathered for the free particle by the Lorentz trans- 
formation itself [7], both the Klein-Gordon and Dirac equations were artificially concocted for the express 
purpose of evading nonlocal integral operators in the configuration representation of relativistic quantum 
mechanics [8, 7]. In light of this fact, it is a truly monumental irony that while the source-free Maxwell 
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equations written in terms of the electric and magnetic fields achieve exactly the configuration-representation 
locality so willfully insisted upon by Klein, Gordon and Dirac, those equations are thereby inconsistent with 
quantization, being far from properly canonically Hamiltonized, and when this issue is duly attended to, 
they utterly lose that "precious" configuration-representation locality to precisely one of the square-root 
operators which Klein, Gordon and Dirac so abjured! An extremely powerful lesson emerges from this: 
there quite simply can be no physically legitimate union of special relativity with quantum theory without 
those square-root operators. The efforts of Klein, Gordon and Dirac to rid relativistic quantum theory of 
these square-root operators and achieve configuration-representation locality merely engenders a dismal list 
of theoretically inappropriate or grotesquely unphysical consequences [7]. The Dirac free particle theory, 
for example, presents a number of operators that are without question physical observables, such as the 
three components of velocity and the energy (i.e., the Dirac Hamiltonian itself), which nonetheless fail to 
mutually commute when the limit % — > is taken! That is grotesquely unphysical, and quite enough to 
permanently consign Dirac theory to the dustbin. But just for good measure, these commutators diverge 
in the nonrclativistic limit c — > oo! Bereft of a Hamiltonian, Klein-Gordon theory has no time evolution 
operator and no Heisenberg picture. For the same underlying reason, it manifests negative probabilities. In 
short, it is so hopelessly crippled that it cannot be regarded as quantum theory at all. These are items from 
the dismal list of Dirac and Klein-Gordon equation shortcomings, but there is not so much as a single item 
in that pejorative list which pertains to the square root Hamiltonianl With this utterly lopsided accounting 
of the theoretical pros and cons, it is far past time for the theoretical physics community to finally awaken 
to just what in its repertoire needs to be revised [7] to give the operator (m 2 c 4 + cp| 2 ) = exactly the same 
standing in relativistic quantum physics that the operator |p| 2 /(2m) properly has in nonrelativistic quantum 
physics. 
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